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ABSTRACT 

Aims. We study the effect of rotation on sheared turbulence due to differential rotation in the solar tachocline. 

Methods. By solving quasi-linear equations for the fluctuating fields, we derive turbulence amplitude and turbulent transport coefficients 
(turbulent viscosity and diffusivity), taking into account the effects of shear and rotation on turbulence. We focus on the regions of the 
tachocline near the equator and the poles where rotation and shear are perpendicular and parallel, respectively. 

Results. For parameter values typical of the tachocline, we show that the shear reduces both turbulence amplitude and transport, more 
strongly in the radial direction (parallel to the shear) than in the horizontal one, resulting in anisotropic turbulence. Rotation further 
reduces turbulence amplitude and transport at the equator whereas it does not have much effect near the pole. The interaction between 
shear and rotation is shown to give rise to a novel non-diffusive flux of angular momentum (known as the A-effect), possibly offering a 
mechanism for the occurrence of a strong shear region in the solar interior. Further implications for the transport in the tachocline are discussed. 
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1. Introduction 

Solar rotation, both global and differential, plays a crucial role 
in the dynamical processes tak ing place in the Sun. For exam- 
ple, m ean-field dynamo theory ( MoffatJl978 : Krause & Ra dler 
Il980h uses these two ingredients to explain the presence of 
the solar magnetic field. Firstly, the differential rotation in the 
tachocline is responsible for the in-situ generation of toroidal 
field via the Q effect. Secondly, convection under the influence 
of (global) rotation leads to helical motion in the convection 
zone which permits the creation of poloidal fields via the a ef- 
fect. 

To understand the differential rotation in the convective 
zone, many authors investigated the structure of turbulence in 
rotating bodies. The main feature of this type of turbulence is 
the appearance of non-diffusive term in the transport of an- 
gular momentum which prevents a solid b ody rotation from 



being a solution o f the Reynolds equation (ILebedinskvll 1 94 1 



Kippenhah^[l963|). Starting from Navier-Stokes equation, it is 
possible to show that these fluxes arise when there is a cause 
of anisotropy in the sys tem, either due to an anisotropic back- 
ground turbulence (see Riidiger 1989, and references therein) 
or due to inhomo geneity such as an underlying stratification 
( Kic hatinovll 19871) . To explain the observed internal differen- 
tial rotation and the depletion of light elements on the surface 
of the Sun, it is also of prime importance to understand the in- 
fluence of rotation on turbulent transport coefficients such as 
the turbulent heat conductivity and the turbulent diffusivity of 
particles. 



The purpose of this Letter is to investigate the effect of 
(global) rotation on the tachocline transport. In our previous 
works, we have studied the turbulent transport in the tachocline 
by taking into account the crucial effect of shearing, the so- 
called sh ear stabilis a tion, due to a strong ra dial differential 
rotation dKiml 120051: ILeprovost & Kiml 120061) . We have also 
incorporated the interaction of this sheared turbulence with 
different types of waves that can be excited in the Sun due 
to magnetic fields ( Leprovost & Kiml 120071) or stratification 
( Kim & Leprovostl2006 ). In this Letter, we elucidate the effect 
of global rotation on sheared turbulence by studying a (local) 
Cartesian model valid near the equator and the poles. We con- 
sider a turbulence driven by an external forcing such as plumes 
from the convection zone and perform a quasi-linear analysis to 
derive the dependence of turbulence amplitude and transport on 
rotation a nd shear. Compared to tw o-dimensional turbulence 
studied in ILeprovost & Kiml d2007l) . global rotation supports 
the propagation of inertial waves in three dimensions (3D), 
which interact with the shear flow, playing an important role in 
the overall turbulent transport. In particular, we report a novel 
result that the momentum transport can be not only due to eddy 
viscosity but also to non-diffusive A-effect in the tachocline. 
Non trivial A-effect here results from an anisotropy induced 
by shear flow on the turbulence even when the driving force is 
isotropic, in contrast to the case without sh ear flow where thi s 
effect exists only for anisotropic turbulence (iKichafinovi 1987). 
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Fig. 1. The configuration in our model: !K, Q and 6 are the 
shearing rate, the rotation rate and the co-latitude, respectively. 

2. Model 

To elucidate the effect of rotation on sheared turbulence, 
we model the tachocline by an incompressible fluid in local 
Cartesian coordinates (x,y,z) in a rotating frame (see Fig. Q]) 
with average rotation rate Q. making an angle with a mean 
shear flow in the azimuthal direction: Uo = -xM J. We study 
the effect of this large-scale shear on the transport properties 
of turbulence by assuming the velocity as a sum of a radial 
shear (i.e. in the x-direction) and fluctuations: u = Uq + v = 
Uo(x)j + v = —xJli + v. We resort to t he quasi-linear approx- 
imation dTownsendll 1 9761: iMoffattll 1978b . also called first-order 
smoothing (FOSA) or second-order correlation approximation 
(SOCA), where the product of fluctuations is neglected to ob- 
tain the following equations for the evolution of the fluctuating 
velocity field: 

d,\ + U () ■ Vv + v • VU = -Vp + vV 2 v + f - £1 x v , (1) 
Vv = 0, 

where p and f are respectively the small-scale components of 
the pressure and forcing, and £1 = 2£2. 

To study the influence of rotation and shear on the parti- 
cle and heat transport, we have to supplement Eq. (Q~|) with an 
advection-diffusion equation for these quantities. We here fo- 
cus on the transport of particles since a similar result also holds 
for the heat transport. The density of particles is also written as 
a sum of a large-scale component Nq and fluctuations n. Using 
the quasi-linear approximation, the fluctuating concentration of 
particles is governed by the following equation: 



d,n + Uo • Vn + v • VjV = DV 2 n , 



(2) 



where D is the molecular diffusivity of particle. Note that, in the 
case of heat equation, D should be replaced by the molecular 
heat conductivity^. 

As the large-scale velocity is in the y (azimuthal) direction, 
we are mostly interested in the momentum transport in that di- 
rection. The equation for the (large-scale) azimuthal velocity 
Uq is then given by Navier-Stokes equation with the contribu- 
tion from fluctuations given by V ■ R, where R = (vv v ) is the 



Reynolds stress. One can formally Taylor expand R with re- 
spect to the gradient of the large-scale flow: 

Ri = AiU () - v T d x U S n + ... = A,f/ + v T tt6 n + ■ . . . (3) 

where we introduced two coefficients A, and vj. The effect of 
the turbulent viscosity vj is simply to change the viscosity from 
the molecular value v to the effective value v + vj. In compari- 
son, the first term A, t/o in Eq. (0) is proportional to the veloc- 
ity rather than its gradient. This means that it does not vanish 
for a constant velocity field and can thus lead to the creation 
of gradient in the velocity fi eld. This term is equivalent to the 



g-effe ct in dynamo theory dParkerlll95 5; Steenbeck & Krause 



1966), and has been known as the A-effect (Krause & Riidigei 



1974) or anisotropic kinetic alpha (AKA)-effect (iFrisch et al 



19871) . Similarly, the transport of species in the large scales is 
governed by an advection diffusion equation where the molec- 
ular diffusivity is supplemented by a turbulent diffusivity, de- 
fined as (v,'«) = -D^djNo. 

To calculate the turbulence amplitude and transport coeffi- 
cients (Reynolds stress and turbulent diffusivity), we prescribe 
the forcing in Eq. (Q~|) to be incompressible, isotropic and short 
correlated in time (modelled by a 5-function) with a power 
spectrum F. Specifically, we assume: 

(/fcki, ?i)/;(k 2 , f 2 )> = Tf (27r) 3 5(k! + k 2 ) S(h - h) x (4) 

FWidij-kikj/e) . 

The angular brackets stand for an average over realisations of 
the forcing, and tj is the (short) correlation time of the forcing. 
We solve Eqs. (HJ and (|2]) in the case of unit Prandtl number 
(v = D) for simplicity, and use the results of these and Eq. (@]l 
to compute the turbulent transport. In this Letter, we focus on 
the case relevant to the Sun and discuss the implications of our 
findings for the dynamics of the tachocline. Note that in the 
tachocline, ~ 3 x 10~ 6 s _1 by using the upper limit on the 
tachocline thickness (5% of the solar radius) while the average 
rotation rate of the interior is Q ~ 2.6 X 10 s" 1 . Therefore, 
the ratio of rotation to shear is very close to unity. However, 
if the tachocline is found to be thinner, the value of the shear- 
ing rate JK becomes larger, with |£2/J?l| < 1. Furthermore, with 



vkyl\Ji\ is a small 



molecular viscosity of v ~ 10 cm s , ^ 
parameter for a broad range of reasonable length scales in the 
azimuthal direction, L y > 10 4 cm. Thus, the present results are 
valid for slow rotation (Q <k |^[|) and strong shear (£ <sc 1). 

3. Near the equator (6 = n/2) 

Near the equator, the rotation and the direction of the shear are 
orthogonal. Expanding the velocity in powers of Q = we 
compute the turbulence amplitude and transport coefficients up 
to second order in |Q| <k 1. 



3.1. Turbulence intensity 

In the stro ng shear limit (g <K 1), using similar algebra as in 
Kim (12005b . we obtain the turbulent intensity in the radial di- 
rection as follows: 



<v?> 



TJL f 

a J 



d 3 k 
(2n) 



■F(k) 



_ k 2 

L (k) + Q^Li(k) 



(5) 
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Here, k^ = k 2 , + k 2 is the (squared) amplitude of the wave num- 
ber in the horizontal plane, and Lq and L\ are two positive defi- 
nite integrals depending only on the wave number k. Note that 
the details of Lq and L\ are not important for our discussions. 
Therefore, the turbulence intensity (v^) in Eq. (O increases for 
£2 > whereas it decreases for £2 < 0. This is because a weak 
rotation destabilises sheared turbulence for fl > whereas it 
stabilises for Q < dBradshawHl969llSalhi & Cambonll 1997th 
For the other components of the turbulence amplitude, we 
obtain: 



<v'> 



2 f 



d 5 k 



F(k)L 2 (k)\- 



1/3 



(6) 



r(l/3) + Q^r(4/3)(-lna 

Here, F is the Gamma function and L 2 is a positive definite 
function of k. Here again, the turbulence amplitude can in- 
crease or decrease depending on the sign of £L Furthermore, the 
correction due to rotation has a logarithmic dependence on the 
shear, contrary to the case of the radial amplitude. Note that the 
anisotropy in the first correction is less pronounced than in the 
leading order term. As a result, the turbulence due to shearing 
becomes less anisotropic. This illustrates the tendency of rota- 
tion to lead to almost isotropic turbulence. An equation similar 
to © is also found for (\ 2 ) with a slightly different function 
L 2 . However, performing the angular integration, we find that 
(v 2 ) is larger than_(v?), in agreement with numerical simula- 
tions ( Lee et al1ll990l) . This effect is present for £2 = and is 
thus caused by the shear only. 

3.2. Transport of angular momentum 

In the strong shear limit (£ <K 1), the transport of angular 
momentum can be found as a sum of two terms (v^v,,) = 
vjJ\+ A A £2, the first and second being even and odd with the ro- 
tation rate, respectively. The first term, the turbulent viscosity, 
takes the following form: 

d 3 k 



Vj 



TJ ~ f 

^ 2 J 



(2n) : 



F(k) 



(7) 



where L3 is a posi tive d e finite function of k. For £2 = we re- 
cover the result of Kim (2005) showing that the turbulent vis- 
cosity is reduced proportionally to J\T 2 for strong shear. The 
turbulent viscosity can be either positive or negative depending 
on the relative magnitude of the two terms inside the integral. 
In the 2D limit (where L3 = 0), we can easily see that the tur- 
bulent viscosity is negative. On the contrary, in the case of an 
isotropic forcing in 3D, the turbulent viscosity is positive. 

The correction due to the rotation is proportional to Q and 
is thus odd in the rotation. This is the so-called A-effect, a non- 
diffusive contribution to Reynolds stress, which can be shown 
to be: 

Ax ~^ J (| f(W)( - h ^ (8) 



It is important to emphasise that this non trivial A-effect re- 
sults from an anisotropy induced by shear flow on the turbu- 
lence even when the driving force is isotropic. This should be 



contrasted to the case without shear flow whe re non-diffusive 
fluxes emerge on ly for anisotropic forcing dRiidigei 1980; 
Kichatinovl fl9861) . In our case, the anisotropy in the velocity 
field is not artificially introduced in the system but is created 
by the shear and calculated self-consistently. 

3.3. Transport of particles 

In the strong shear limit (£ <k 1), the computation of (nv) gives: 

k 2 



Di 



IL f^l 
& J (2tt 

T f r (fl 

~& J (2^ 



F(k)L 5 (k) 



1 + £2 



■ln£ 



-F(k)L 6 (k)\ — 



k 2 3 



■ln£ 



(9) 



,(10) 



where L5 and are positive definite functions. Eqs. (l9lfT0li 
show that the effect of rotation on the transport of particles de- 
pends on the sign of £2: for £2 > 0, the transport is increased 
whereas it is reduced when £2 < 0. This is again because a weak 
rotation destabilises sheared turbulence for £2 > whereas it 
stabilises for £2 < 0. Note that a similar behaviour was also 
found in turbulence intensity, given in Eqs. © and ©. The 
correction term due to rotation in D" and Dj in Eq. d9lfT0b de- 
pends weakly on the shear by a logarithmic factor | ln£|, which 
cannot be too large even for ^ « 1, and is of the same order for 
transport in different directions. Thus, the scaling of the turbu- 
lent diffusivity is roughly the same as that in the case without 
rotation: the radial transport (oc JT 2 ) is more reduced than the 
horizontal one (oc j?T 4 / 3 ). This result should be contrasted to 
the rapid rotation limit where the transport in the radial direc- 
tion was larger (bu t only by a fac tor 2) than the one in the hor- 
izontal direction (Rudiger 19891) . These results thus highlight 
the crucial role of shear in transport, in particular, by introduc- 
ing anisotropy. 

4. Near the poles (6 = 0) 

Near the poles, the directions of rotation and shear can be taken 
to be parallel. Contrary to the case at the equator, we find that 
the leading order correction (proportional to £2) vanishes in the 
case of an isotropic forcing (due to the fact that these terms 
are odd in one of the components of the wave number) for all 
the previously calculated quantities. Consequently, the turbu- 
lent amplitude, viscosity and diffusivity are given by Eqs. ([5} 
UK © and d9l[TCfl i with £1 = 0, respectively. Similarly, the 
Reynolds stress involving the radial component ((v x v y )) van- 
ishes. However the component of the Reynolds stress involv- 
ing the latitudinal velocity ((VyV,)) does not vanish and is odd 
in Q. Thus, the A-effect appears here in (v y \ z ) (recall that at the 
equator, the A-effect was present only in (V jV^,)), and takes the 
following form: 



A- 



T f r d 3 ) 

(27 



2 k I 3 \ 2/3 
f m \2i) [i7(k) ' 



(id 



Eq. ( fTTT i shows that A z is of indefinite sign, as L 7 and L% are 
two positive definite functions but appear with different signs. 
However, as for Eq. ((8), in the case of an isotropic forcing, the 
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function L 7 dominates over leading to a negative A z . This 
A-effect arises even in the case of an isotropic forcing as the 
shear favours fluctuations in the y-direction compared to that in 
the z-direction (see S3. Il l, leading to anisotropic turbulence in 
the plane perpendicular to the rotation. Eq. ( fTTb also shows that 
A z is larger (scaling as JT 4 ' 3 ) than A x in the equatorial case. 

5. Implications for the tachocline 

Table Q] summarises the scaling of our results with the shear 
y\. and rotation Q. in the limit of weak rotation (O <k \^l\) and 
strong shear (£ = v&f\$l\ «C 1). 





Equator 


Pole 


<v?> 


ft-' h+cCl] 


ft- 1 


<V?> ~ < V 2> 




1 +cn|lnf|] 


ft- 2 ' 3 


y T 


ft- 2 


ft- 2 


A x 


ftT 2 \ki£\ 





A; 





-ft- 4 ' 3 


Df 


ft- 2 


l + cn|lnf|] 


ft- 2 


~ Li- 


ft- 4 ' 3 


1 +cn|lnf|] 


ft- 4 ' 3 



Table 1. Scaling of turbulence amplitude, turbulent viscosity 
(vt), A-effect and turbulent diffusivity with the shear J{ and 
rotation Q. — Q/^[. c is a positive constant of order unity. 



As both vj and A are positive, J{ must be negative. 
Interestingly, this is exactly what the observation indicates (a 
rotation rate decreasing towards the interior at the equator). We 
cannot use such a simple equation to predict the sign of the 
shear at the pole as the A-effect now appears in the (v v v,) com- 
ponent of the Reynolds stress and is found to scale as , 
which is larger than A A (oc J?T 2 ) near the equator. 

Furthermore, we find that the transport of chemical species 
is reduced by shear and more severely in the the radial direc- 
tion than the horizontal one (by a factor J?T 2 and J{~ 4 ' 3 respec- 
tively), with a further, but weak reduction due to rotation at 
the equator only. As a result, the radial mixing (in the x direc- 
tion) of chemicals is slightly more efficient at the pole than at 
the equator. These results thus suggest that turbulent mixing of 
light elements is weak (due to shearing) in the tachocline and 
also that the depletion of light elements on the solar tachocline 
depends on latitude. However, this latitudinal dependency is 
unlikely to be observed on the surfa ce due to the rapid mixing 
in the convection zone dSpruitlll977h . Since a similar result ap- 
plies to heat transport, a faster heat transport at the pole could 
contribute to making hotter poles than equator at the solar sur- 
face. Finally, we note that all our results he re are valid in the 
absen ce of stratification and magnetic fields dKim & Leprovostl 
2006) and will extend our study to include these in future pub- 
lications. Of particular interest would be the sign of vj and A, 
and thus the resulting sign of M. 
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The first two rows show that, both near the equator and 
poles, the turbulence amplitude in the radial direction is more 
reduced by the shear than in the horizontal one, by a factor of 
^T 1 and J{- 2 ' 3 , respectively. These results thus imply an ef- 
fectively stronger turbulence in the horizontal (y-z) plane than 
the one in the radial (x) direction. Furthermore, at the equator, 
Q < as the shear increases from the interior towards the con- 
vective zone. Thus, the turbulence amplitude is further reduced 
by the rotation in all directions, but the turbulence in the ra- 
dial direction is more suppressed than the one in the horizontal 
direction by a factor ln£. 

A novel result is also found for the Reynolds stress, which 
involves two contributions: the turbulent viscosity and the A- 
effect. The latter, a source of non-diffusive flux, is present 
even with isotropic forcing due to the shear-induced anisotropy 
dKimll2005l) . Note that for an isotropic forcing, the turbulent 
viscosity is found to be positive and therefore cannot act on 
its own as a source of differential rotation in the tachocline. 
The radial differential rotation can, however, be created by the 
A-effect: since it does not depend only on the gradient of an- 
gular velocity, it does not vanish for a uniform rotation and 
thus prevents the uniform rotation from being a solution of the 
large-scale momentum equation. To determine the sign of 
as a result of the A-effect, we seek a solution of the large-scale 
turb ulent equations by demanding the Reynolds stress to vanish 
(e.g. Kic hatinovll 19871) . In the equatorial case, we then obtain 
the following equation for the shear: 

vrcsrm + Ac?on = o . (12) 



References 

Bradshaw, P. 1969, J. Fluid Mech., 36, 177 
Frisch, U., She, Z. S., & Sulem, P. L. 1987, Physica D, 28, 382 
Kichatinov, L. L. 1986, Geophys. Astrophys. Fluid Dyn., 35, 
93 

Kichatinov, L. L. 1987, Geophys. Astrophys. Fluid Dyn., 38, 

273 

Kim, E. 2005, A&A, 441, 763 

Kim, E. & Leprovost, N. 2006, submitted to A&A 

Kippenhahn, R. 1963, ApJ, 137, 664 

Ki-ause, F. & Radler, K.-H. 1980, Mean field MHD and dynamo 

theory (Pergamon press) 
Krause, F. & Rudiger, G. 1974, Astron. Nachr., 295, 93 
Lebedinsky, A. I. 1941, AZh, 18, 10 

Lee, J. M., Kim, J., & Moin, P. 1990, J. Fluid Mech., 216, 561 

Leprovost, N. & Kim, E. 2006, A&A, 456, 617 

Leprovost, N. & Kim, E. 2007, ApJ, 654, 1 166 

Moffatt, H. K. 1978, Magnetic field generation in fluids (CUP) 

Parker, E. N. 1955, ApJ, 122, 293 

Rudiger, G. 1980, Geophys. Astrophys. Fluid Dyn., 16, 239 
Rudiger, G. 1989, Differential rotation and stellar convection 

(Gordon and Breach) 
Salhi, A. & Cambon, C. 1997, J. Fluid Mech., 347, 171 
Spruit, H. C. 1977, A&A, 55, 151 

Steenbeck, M. & Krause, F. 1966, Z. Naturforsch., Teil A, 21, 
1285 

Townsend, A. A. 1976, The structure of turbulent shear flow, 
2nd edn. (CUP) 



